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Abstract. Compelling evidence for the creation of a new form of matter has been 
claimed to be found in Pb+Pb collisions at SPS. We discuss the uniqueness of often 
proposed experimental signatures for quark matter formation in relativistic heavy ion 
collisions. It is demonstrated that so far none of the proposed signals like J/\\t me- 
son production/suppression, strangeness enhancement, dileptons, and directed flow 
unambigiously show that a phase of deconfmed matter has been formed in SPS Pb+Pb 
collisions. We emphasize the need for systematic future measurements to search for 
simultaneous irregularities in the excitation functions of several observables in order 
to come close to pinning the properties of hot, dense QCD matter from data. 



1. Introduction 



In the last few years researchers at Brookhaven and CERN have succeeded to measure a 
wide spectrum of observables with heavy ion beams, Au+Au and Pb + Pb. While these 
programs continue to measure with greater precision the beam energy-, nuclear size-, and 
centrality dependence of those observables, it is important to recognize the major mile- 
stones passed thusfar in that work. The experiments have conclusively demonstrated the 
existence of strong nuclear A dependence of, among others, J/\\t and \|/' meson production 
and suppression, strangeness enhancement, hadronic resonance production, stopping and 
directed collective transverse and longitudinal flow of baryons and mesons - in and out of 
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the impact plane, both at AGS and SPS energies -, and dilepton-enhancement below and 
above the p meson mass. These observations support that a novel form of "resonance mat- 
ter" at high energy- and baryon density has been created in nuclear collisions. The global 
multiplicity and transverse energy measurements prove that substantially more entropy is 
produced in A + A collisions at the SPS than simple superposition of A x pp would imply. 
Multiple initial and final state interactions play a critical role in all observables. The high 
midrapidity baryon density (stopping) and the observed collective transverse and directed 
flow patterns constitute one of the strongest evidence for the existence of an extended pe- 
riod (At w 10 fm/c) of high pressure and strong final state interactions. The enhanced 
\|/ suppression in S + U relative to p + A also attests to this fact. The anomalous low 
mass dilepton enhancement shows that substantial in-medium modifications of multiple 
collision dynamics exist, probably related to in-medium collisional broadening of vector 
mesons. The non-saturation of the strangeness (and anti-strangeness) production shows 
that novel non-equilibrium production processes arise in these reactions. Finally, the cen- 
trality dependence of J/\\t absorption in Pb + Pb collisions presents further hints towards 
the nonequilibrium nature of such reactions. Is there evidence for the long sought-after 
quark-gluon plasma that thusfar has only existed as a binary array of predictions inside 
teraflop computers? 

As we will discuss, it is too early to tell. Theoretically there are still too many "sce- 
narios" and idealizations to provide a satisfactory answer. Recent results from microscopic 
transport models as well as macroscopic hydrodynamical calculations differ significantly 
from predictions of simple thermal models, e. g. in the flow pattern. Still, these nonequilib- 
rium models provide reasonable predictions for the experimental data. We may therefore 
be forced to rethink our concept of what constitutes the deconfined phase in ultrarelativis- 
tic heavy-ion collisions. Most probably it is not a blob of thermalized quarks and gluons. 
Hence, a quark-gluon plasma can only be the source of differences to the predictions of 
these models for hadron ratios, the J/\\t meson production, dilepton yields, or the exci- 
tation function of transverse flow. And there are experimental gaps such as the lack of 
intermediate mass A « 100 data and the limited number of beam energies studied thus- 
far, in particular between the AGS and SPS. Now the field is at the doorstep of the next 
milestone: A + A at y/s — 30 — 200 AGeV which have started a few months ago. 

2. JA|/ production 

The QCD factorization theorem is used to evaluate the PQCD cross sections of heavy 
quarkonium interactions with ordinary hadrons. However, the charmonium states (here de- 
noted X) are not sufficiently small to ignore nonperturbative QCD physics. Thus, we eval- 
uate the nonperturbative QCD contribution to the cross sections of charmonium-nucleon 
interaction by using an interpolation between known cross sections [ |[]. The J/\\t-N cross 
section evaluated in this paper is in reasonable agreement with SLAC data [ Eft. 

Indeed, the A-dependence of the J/\\t production studied at SLAC at Ej nc ~ 20 GeV 
exhibits a significant absorption effect [ Q leading to a a b s (J /\\i-N) = 3.5 ±0.8 mb. It 
was demonstrated [ ||] that, in the kinematic region at SLAC, the color coherence effects 
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are still small on the internucleon scale for the formation of //\|/'s. So, in contrast to 
the findings at higher energies, at intermediate energies this process measures the genuine 
J/\\f-N interaction cross section at energies of ~ 15-20 GeV [ M. 

To evaluate the nonperturbative QCD contribution we use an interpolation formula [ 
[|] for the dependence of the cross section on the transverse size b of a quark-gluon config- 
uration Three reference points are used to fix our parametrization of the cross sections (cf. 
Tab. [I]). The X-N cross sections is calculated via: o = / o(b) ■ |*P(;t,)>,z)| 2 cbcd)>dz, where 
*P(x,y,z) is the charmonium wave function. In our calculations we use the wave functions 
from a non-relativistic charmonium model (see [ ^]). 

Table 1. The total quarkonium-nucleon cross sections a. For the % two values arise, due to the spin 



cc-state 


JA|/ 


¥' 


XclO 


XcU 


O (mb) 


3.62 


20.0 


6.82 


15.9 



We follow the analysis of [ [7j to evaluate the fraction of J/\\f's (in pp collisions) that 
come from the decays of the % and \|/'. So, the suppression factor S of 7/\|/'s produced in 
the nuclear medium is calculated as: 

5 = 0.6- (0.92 • S J /V + 0.08 -S^)+ 0.4 • Here S x are the respective suppression factors 
of the different pure charmonium states X in nuclear matter. The S x are for minimum bias 
pA collisions within the semiclassical approximation (cf. [ 

The charmonium states are produced as small configurations, then they evolve to their 
full size. Therefore, if the formation length of the charmonium states, If, becomes larger 
than the average internucleon distance, one has to take into account the evolution of the 
cross sections with the distance from the production point [ || . 

The formation length of the is given by If w 2/?/(m 2 , — rrij^), where p is the 
momentum of the J/\\t in the rest frame of the target. For a J/\\t produced at midrapidity at 
SPS energies, this yields If w 3 fm. Due to the lack of better knowledge, we use the same 
If « 3 fm for the %. For the \\t' we use If w 6 fm, because it is not a small object, but has 
the size of a normal hadron, i.e. the pion. For E[ a b = 800 AGeV we get a factor of two for 
the formation lengths due to the larger Lorentz factor. 

However, this has a large impact on the \|/' to 7/\|/-ratio depicted in Fig. [j], which shows 
the ratio 0.019 -Sy'/Sr/y calculated with (squares (200 GeV) and triangles (800 GeV)) and 
without (crosses) expansion. The factor 0.019 is the measured value in pp collisions, be- 
cause the experiments do not measure the calculated value Sy /Sj /y but (_B w ,o(\|/ ) ) / (B^a ) . 
B w , are the branching ratios for J/y, \|/' — > /jp. 

The calculations which take into account the expansion of small wave packages show 
better agreement with the data (circles) (taken from [ ||]) than the calculation without ex- 
pansion time, i.e. with immediate J/\\t formation, If = 0. We calculated this effect both at 
Eiab = 200AGeV and 800 AGeV. The data have been measured at different energies (£/„/, 
= 200, 300, 400, 450, 800 GeV and ^ = 63GeV). One can see that this ratio is nearly 
constant in the kinematical region of the data, but it decreases at smaller momentum (e.g. 
E[ u b = 200 AGeV and y < 0) due to the larger cross section of the \|/. 

However, the P-states yield two vastly different cross sections (see Tab. [l]) for Xio and 
%n, respectively. This leads to a higher absorption rate of the %u as compared to the 
Xio- This new form of color filtering is predicted also for the corresponding states of other 
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Fig. 1. Left: The ratio 0.019 ■ Sy /S// v is shown in /?A (crosses) in comparison to the data 

(circles). The squares and the triangles shows the ratio calculated with the expansion of 

small wave packages. Right: The ratio of J/\\i to Drell-Yan production as a function of Et 

for Pb+Pb at 160 GeV. . . . , •,«.,/, 

hadrons; e.g. for the bottomium states which are proposed as contrast signals to the J/\y s 

atRHICandLHC! 

Furthermore it is important to also take into account comoving mesons. Therefore 
we use the UrQMD model [ [k^, 11 ]. Particles produced by string fragmentation are not 
allowed to interact with other hadrons - in particular with a charmonium state - within their 
formation time (on average, Tjr w 1 fm/c). However, leading hadrons are allowed to interact 
with a reduced cross section even within their formation time . The reduction factor is 1/2 
for mesons which contain a leading constituent quark from an incident nucleon and 2/3 for 
baryons which contain a leading diquark. 

Figure [j] shows the J/\\i to Drell-Yan ratio as a function of Et for Pb+Pb interactions at 
160 GeV compared to the NA50 data [||, |||]. The normalization of B^a(J /\\f) /c(DY) 
46 in pp interactions at 200 GeV has been fit to S+U data within a geometrical model [ 

The application of this value to our analysis is not arbitrary: the model of Ref. [ 
renders the identical Et -integrated J/\\t survival probability, S = 0.49, as the UrQMD cal 
culation for this system. An additional factor of 1 .25 [ [l4]] has been applied to the Pb+Pb 
calculation in order to account for the lower energy, 160 GeV, since the J/\\t and Drell-Yan 
cross sections have different energy and isospin dependencies. 

The gross features of the Et dependence of the J/\\t to Drell-Yan ratio are reasonably 
well described by the model calculation. No discontinuities in the shape of the ratio as 
a function of Et are predicted by the simulation. The new high Et data [ [j~3|] decreases 
stronger than the calculation. This could be caused by underestimated fluctuations of the 
multiplicity of secondaries in the UrQMD model. This occurs, since high £Y-values are a 
trigger for very central events with a secondary multiplicity larger than in average [ Oh. 



I 



3. Dilepton production 



Beside results from hadronic probes, electromagnetic radiation - and in particular dilep- 
tons - offer an unique probe from the hot and dense reaction zone: here, hadronic matter is 
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almost transparent. The observed enhancement of the dilepton yield at intermediate invari- 
ant masses (M e + e - > 0.3 GeV) received great interest: it was prematurely thought that the 
lowering of vector meson masses is required by chiral symmetry restoration (see e.g.[ |l6|| 
for a review). 



Fig. 2. 

Microscopic calculation of the dilepton pro- 
duction in the kinematic acceptance region of 
the CERES detector for Pb+Au collisions at 
158 GeV. No in-medium effects are taken into 
account. Plotted data points are taken at CERES 
in '95. 




Fig. H shows a microscopic UrQMD calculation of the dilepton production in the kine- 
matic acceptance region of the CERES detector for Pb +Au collisions at 158 GeV. This is 
compared with the '95 CERES data[ |||. Aside from the difference at M sa 0.4 GeV there 
is a strong enhancement at higher invariant masses. It is expected that this discrepancy at 
m > 1 GeV could be filled up by direct dilepton production in meson-meson collisions [ [19| 
as well as by the mechanism of secondary Drell-Yan pair production proposed in[ 



4. Strangeness production 

Strange particle yields are most interesting and useful probes to examine excited nuclear 
matter [ ^jl ^2|, ^4], 26 , E7|] and to detect the transition of (confined) hadronic 
matter to quark-gluon-matter. The relative enhancement of strange and especially multi- 
strange particles in central heavy ion collisions with respect to peripheral or proton induced 
interactions have been suggested as a signature for the transient existence of a QGP-phase 
[ |TJ|. Here the main idea is that the strange (and antistrange) quarks are thought to be pro- 
duced more easily and hence also more abundantly in such a deconfined state as compared 
to the production via highly threshold suppressed inelastic hadronic collisions. The rela- 
tive enhanement of (anti)hyperons has clearly been measured by the WA97 an the NA49 
collaboration in Pb-Pb collisions as compared to p-Pb collisions [ ^5[ ^6|] . This data has 
been investigated within microscopic transport models (e.g. UrQMD [ |10|]). In [ |5J Q 
it was shown that within such an appraoch strangeness enhancement is predicted for Pb- 
Pb due to rescattering. However, for central Pb-Pb collisions the experimentally observed 
hyperon yields are underestimated by the calculation in [ |2[ E8I]. This result seems to 
confirm the conclusion that a deconfined QGP is formed in Pb-Pb collisions at SPS. But 
in [ ^7|, [29]] it was shown, that the antihyperon production by multi-mesonic reactions like 
n\K + 112K — ¥ Y +p could drive these rare particles towards local chemical quilibrium with 
pions, nucleons and kaons on a timescale of 1-3 fm/c. Accordingly this mechanism, which 
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is a consequence of detailed balance could provide a convenient explanation for the antihy- 
peron yields at CERN-SPS energies without any need of a deconfined quark-gluon-plasma 
phase. At the moment such back-reactions cannot be handled within the present transport 
codes. Therefore the aim for the future will be to find a way to include these processes in 
microscopic transport models. 



5. Particle ratios 

Ideal gas model calculations have been used for a long time to calculate particle produc- 
tion in relativistic heavy ion collisions, see e.g. [ 5l[ ^2], |33|, |34], |3^j. Fitting the 
particle ratios as obtained from those ideal gas calculations to the experimental measured 
ratios at SIS, AGS and SPS for different energies and different colliding systems yields a 
curve of chemical freeze-out in the T — /j plane. Now the question arises, how much the 
deduced temperature and chemical potentials depend on the model employed. Especially 
the influence of changing hadron masses and effective potentials should be investigated, 
as has been done for example in [ |3(], |7l [38[ |3^] . This is of special importance for the 
quest of a signal of the formation of a deconfined phase, i.e. the quark-gluon plasma. As 
deduced from lattice data [ the critical temperature for the onset of a deconfined phase 
coincides with that of a chirally restored phase. Chiral effective models of QCD therefore 
can be utilized to give important insights on signals from a quark-gluon plasma formed in 
heavy-ion collisions. 

Therefore we compare experimental measurements for Pb+Pb collisions at SPS with 
the ideal gas calculations and results obtained from a chiral SU(3) model [ g 0]. This 
effective hadronic model predicts a chiral phase transition at T w 150Me V. Furthermore the 
model predicts changing hadronic masses and effective chemical potentials, due to strong 
scalar and vector fields in hot and dense hadronic matter, which are constrained by chiral 
symmetry from the QCD Lagrangean. 



In [ 32 1 the ideal gas model was fitted to particle ratios measured in Pb+Pb collisions 
at SPS. The lowest % 2 is obtained for T = 168MeV and fi q = 88.67MeV. Using these 
values as input for the chiral model leads to dramatic changes due to the changing hadronic 
masses in hot and dense matter [ ^] and therefore the freeze-out temperature and chemical 
potential have to be readjusted to account for the in-medium effects of the hadrons in the 
chiral model. We call the best fit the parameter set that gives a minimum in the value 

of x 2 , with % 2 = Y,i — ^ • Here rf xp is the experimental ratio, r" mdel is the ratio 

calculated in the model and a,- represents the error in the experimental data points as quoted 
in [ B2p. The resulting values of % 2 for different T — /u pairs are shown in figure [3| In 
all calculations /j s was chosen such that the overall net strangeness f s is zero. The best 
values for the parameters are T = 144MeV and /j c/ sw 95MeV. While the value of the 
chemical potential does not change much compared to the ideal gas calculation, the value 
of the temperature is lowered by more than 20 Me V. Furthermore Figure || shows, that the 
dropping effective masses and the reduction of the effective chemical potential make the 
reproduction of experimentally measured particle ratios as seen at CERN's SPS within this 
model impossible for T > T c . Using the best fit parameters a reasonable description of the 
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O SPSPb+Pb 

taken from Braun-Munzinger et al, PLB 465(1999) 15 
A chiral T=144 McV, /j q =95 McV 
□ ideal gas calculation, PLB 465(1999) 15 




T [Me V| 



Fig. 3. X 2 (left) an d resulting particle ratios compared to ideal gas calculation and data (rigft) for 



chiral model, data taken from [|32p. The best fit parameters are T = 144MeV and fj q ~ 95MeV. 

particle ratios used in the fit procedure can be obtained (see fig.|3|, data from [ f32|]). 

We want to emphasize, that in spite of the strong assumption of thermal and chemical 
equilibrium the obtained values for T and /j differ significantly depending on the under- 
lying model, i.e. whether and how effective masses and effective chemical potentials are 
accounted for. Note that we assume implicitly, that the particle ratios are determined by the 
medium effects and freeze out during the late stage expansion - no flavor changing colli- 
sions occur anymore, but the hadrons can take the necessary energy to get onto their mass 
shall by drawing energy from the fields. Rescattering effects will alter our conclusion but 
are presumably small when the chemical potentials are frozen. 



6. Collective flow and the EOS 

The in-plane flow has been proposed as a measure of the "softening" of the EoS [ [|3|], there- 
fore we investigate the excitation function of directed in-plane flow. A three-fluid model 
with dynamical unification of kinetically equilibrated fluid elements is applied [ This 
model assumes that a projectile- and a target fluid interpenetrate upon impact of the two 
nuclei, creating a third fluid via new source terms in the continuity equations for energy- 
and momentum flux. Those source terms are taken from energy- and rapidity loss mea- 
surements in high energy pp-collisions. The equation of state (EoS) of this model assumes 
equilibrium only in each fluid separately and allows for a first order phase transition to a 
quark gluon plasma in fluid 1, 2 or 3, if the energy density in the fluid under consideration 
exceeds the critical value for two phase coexistence. Pure QGP can also be formed in every 
fluid separately, if the energy density in that fluid exceeds the maximum energy density for 
the mixed phase. Integrating up the collective momentum in x-direction at given rapidity, 
and dividing by the net baryon number in that rapidity bin, we obtain the so-called directed 
in-plane flow per nucleon. 

Its excitation function (Fig. |j) shows a local minimum at 8 AGeV and rises until a 
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maximum around 40 AGeV is reached. Fig. |] shows the excitation function of directed 
flow calculated in the three-fluid model in comparison to that obtained in a one-fluid cal- 
culation. Due to non-equilibrium effects in the early stage of the reaction, which delay the 
build-up of transverse pressure [pq ], the flow shifts to higher bombarding energies. While 
measurements of flow at AGS[ |46| ] have found a decrease of directed flow with increasing 
bombarding energy, a minimum has so far not been observed. 

In a recent investigation of the directed flow excitation functions [ ^7|] it has been 
shown, that the directed flow excitation functions are sensitive to the underlying EoS and 
that a different EoS can predict a slowly and smoothly decrease of the averaged directed 
flow as a function of bombarding energies. This different behaviour is due to the different 
phase transitions in the underlying equations of state. While in the two phase EoS based 
on a a — co model for the hadronic phase and a bag model for the deconfined phase a first- 
order phase transition occurs, the EoS in [ ^7[| provides a continues phase transition of the 
cross-over type. 



x 1 fluid, no PT 

O 1 fluid, wilh PT 

O 3 fluid, with PT 

A EOS 





IWGeV] [AGeV] 

Fig. 4. (Left) Excitation function of transverse flow as obtained from three fluid hydro- 
dynamics with a first order phase transition and (Right) the slope of the directed in-plane 

momentum per nucleon at midrapidity. 

The slope of the directed in-plane momentum per nucleon at midrapidity, d(p x /N)/dy, 

is shown in Fig. [| as a function of beam energy. We find a steady decrease of d(p x /N) / dy 
up to about top BNL-AGS energy, where the flow around midrapidity even becomes nega- 
tive due to preferred expansion towards p x ■ pi on g < 0. At higher energy, Euib — 40A GeV, 
the isentropic speed of sound becomes small and we encounter the following expansion pat- 
tern : flow towards p x ■ pi ong < can not build up ! Consequently, d(p x /N)/dy increases 
rapidly towards Eu,b = 20 — 40A GeV, decreasing again at even higher energy because of 
the more forward-backward peaked kinematics which is unfavorable for directed flow. 

Thus, the Pb + Pb collisions (40 GeV) runs performed recently at the CERN-SPS may 
provide a crucial test of the picture of a quasi-adiabatic first-order hadronization phase 
transition at small isentropic velocity of sound. 



7. Collective Flow at RHIC 

Let us now compare the first results on elliptic flow (V2) at ^/sjvjv = 130GeV as reported by 
the STAR-Collaboration [ Q with a string hadronic model simulation: The experimental 
data indicates a strongly rising V2 as a function of p t with an average V2 value of 6% at 
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midrapidity and p t approximately 600MeV. While the strong increase of v% with p t has 
been predicted by the UrQMD model [ E9] the absolute magnitude of V2 at p t = 600MeV 
is underpredicted by a factor 3 (cf. fig. 
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Au+Au, s 1/2 — 200 AGeV, min. bias, cut: y cm +l 



Fig. 5. 

Elliptic flow parameter v% at midrapidity as a 
function of transverse momentum in minimum 
biased Au+Au reactions at \fs = 200AGeV 



).0 0.2 0.4 0.6 0.8 1.0 1.2 
Pt (GeV/c) 

When the formation time of hadrons in the initial strings is strongly reduced (to mimic 
short mean free paths in the early interaction region) the calculated flow values approach the 
hydrodynamic limit [ [5(], and get in line with the measured elliptic flow values. This 
shows, that the pressure in the reaction zone is much higher than expected from simple 
stringlike models and supports the breakdown of pure string hadronic dynamics in the 
initial stage of Au-Au-collisions at RHIC energies. However, to get a consistent picture 
and to finally rule out the string hadronic approach the vi values and transverse momentum 
spectra [ |T]] as given by the model calculation need to be exceeded by the experimental 
data. 



8. Insights from quark molecular dynamics 

Further insights about the possible formation of deconfined matter can be obtained from the 
Quark Molecular Dynamics Model (qMD) [ [52]] which explicitly includes quark degrees 
of freedom. The qMD can provide us with detailed information about the dynamics of the 
quark system and the parton-hadron conversion. Correlations between the quarks clustering 
to build new hadrons can be studied [p3fa. 

Figure || shows (for S+Au collisions at SPS energies of 200GeV/N) the number dis- 
tribution for the mean path travelled by quarks forming a hadron (a) from the same initial 
hadron (solid line) and (b) from different initial hadrons (dotted line). 

A measure of the relative mixing within the quark system and thus for thermalization is 
the relative number of hadrons formed by quarks from the same initial hadron correlation 
versus hadrons formed by quarks from different initial hadron correlations. This ratio is 
r = 0.574 ± 0.008 for the S+Au collision. Since a value of r = 1 would indicate complete 
rearrangement of quarks and thus complete loss of correlations in the quark system, one 
would expect a much larger value of r, considering the presumed transition to the quark- 
gluon plasma in Pb+Pb collisions at 160GeV/N, 
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Fig. 6. 

Hadronization in S+Au collisions at SPS (200 
GeV/N): Number density distribution of mean 
diffusion path of quarks forming a hadron from 
the same initial hadron (solid line) and from 
different initial hadrons (dashed line) within 
qMD. Fitting the decay profiles yields diffusion 
lengths of 2.2fm and 4.8fm, respectively. 



The latest data of CERN/SPS on flow, electro-magnetic probes, strange particle yields 
(most importantly multistrange (anti-)hyperons) and heavy quarkonia will be interesting 
to follow closely. Simple energy densities estimated from rapidity distributions and tem- 
peratures extracted from particle spectra indicate that initial conditions could be near or 
just above the domain of deconfinement and chiral symmetry restoration. Still the quest for 
an unambiguous signature remains open. 

Directed flow has been discovered - now a flow excitation function, filling the gap 
between 10 AGeV (AGS) and 160 AGeV (SPS), will be extremely interesting: look for 
the softening of the QCD equation of state in the coexistence region. The investigation of 
the physics of high baryon density (e.g. partial restoration of chiral symmetry via prop- 
erties of vector mesons) has been pushed forward by the 40 GeV run at SPS. Also the 
excitation function of particle yield ratios (jt/p,d/p,K/n...) and, in particular, multistrange 
(anti-)hyperon yields, can be a sensitive probe of physics changes in the EoS. The search 
for novel, unexpected forms of matter, e.g. hypermatter, strangelets or even charmlets 
is intriguing. Such exotic QCD multi-meson and multi-baryon configurations would ex- 
tend the present periodic table of elements into hitherto unexplored dimensions. A strong 
experimental effort should continue in that direction. 

Now we have entered the exciting RHIC era, where the predicted deconfined and chi- 
rally restored phase should be formed and live long enough to produce clear and unambi- 
gious signals of it's existence. The LHC-program will top this scientific endeavour in 4 
years. 
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